Europhysics Letters 



PREPRINT 



Small scale statistics of viscoelastic turbulence 

S. Berti^, a. Bistagnino^ G. Boffetta^ a. Celani^ and S. Musacchio^ 

^ Dipartimento di Fisica Generale and INFN, Universitd degli Studi di Torino, Via 
Pietro Giuria 1, 10125, Torino, Italy 
and ISAG-GNR, Sezione di Torino, Italy 

^ INLN-GNRS, 1361 route des Lucioles, Sophia Antipolis, 06560 Valbonne, France. 



PACS. 47.27.Gs - Isotropic turbulence; homogeneous turbulence. 

PACS. 47.27.E Turbulence simulation and modeling. 

PACS. 47.57.Ng - Polymers and polymer solutions. 



Abstract. - The small scale statistics of homogeneous isotropic turbulence of dilute polymer 
solutions is investigated by means of direct numerical simulations of a simplified viscoelastic 
fluid model. It is found that polymers only partially suppress the turbulent cascade below the 
Lumley scale, leaving a remnant energy flux even for large elasticity. As a consequence, fluid 
acceleration in viscoelastic flows is reduced with respect to Newtonian turbulence, whereas its 
rescaled probability density is left unchanged. At large scales the velocity field is found to be 
unaffected by the presence of polymers. 



The addition of small amounts of long chain polymers produces dramatic effects on flow- 
ing fluids, the most renowned being the reduction of friction drag in turbulent flow at high 
Reynolds numbers [1,2]. Most studies focused on dilute polymer solutions in channel or pipe 
geometry, where boundary effects are important [3], but recent experimental [4,5] and nu- 
merical [6-8] works have shown that polymers affect the turbulent flow even far from (or 
in absence of) boundaries. In particular, Ref. [7] studied the modification of the turbulent 
cascade induced by polymers in numerical simulations of homogeneous, isotropic turbulence. 
In this paper we investigate the effects of polymer addition to the small scale statistics in 
fully developed homogeneous-isotropic turbulence by means of direct numerical simulations 
of a simplified viscoelastic model. We show that, by increasing the elasticity of polymers, 
the energy flux in the turbulent cascade is partially suppressed and transferred to the elastic 
degrees of freedom. This suppression remains partial even for large values of elasticity: as a 
consequence the energy flux to small scales remains finite and the small scale statistics, such 
as acceleration probability density function (pdf), retain some characteristics of Newtonian 
flows. 

The mechanism by which dilute polymer solutions can influence turbulent flows is the ex- 
treme extensibility of polymers. Polymers, typically composed by a large number of monomers, 
at equilibrium are coiled in a ball of radius Rq. In presence of a nonhomogeneous flow, the 
molecule is deformed in an elongated structure characterized by its end-to-end distance R 
which can be significantly larger than Rq. The deformation of molecules is the result of the 
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competition between the stretching induced by differences of velocities and the entropic re- 
laxation of polymers to their equilibrium configuration. Experiments with DNA molecules [9] 
show that this relaxation is linear, provided that the elongation is small compared with the 
maximal extension R <^ Rmax, and can be characterized by a typical relaxation time r [10]. 

These ingredients lead to the simplest model which describes the behavior of a polymer 
in a flow, the dumbbell model. Since in applications the typical size of polymers is smaller 
than the viscous scale of turbulence, stretching is due to velocity gradients and the end-to-end 
distance evolves according to: 



-^iyurR--R+sj-^^ (1) 

where ^ is a Brownian process with correlation {^i{t)^j{t')) = 5ij5{t — t'). 

The relative importance between polymer relaxation and stretching is measured by the 
Weissenberg number VFi, defined as the product of r and the characteristic velocity gradient. 
When Wi <^ 1 relaxation is fast compared to the stretching time and polymers remain in 
the coiled state. For Wi ^ 1, on the contrary, polymers are substantially elongated. The 
transition point is called the coil-stretch transition and occurs at Wi = 0(1). 

In the case of dilute solutions, for which the polymer concentration n satisfies uRq 1, 
the influence of polymers in the coiled state on the fluid is negligible. Above the coil-stretch 
transition, polymers start to affect the flow. This regime is characterized by large elongations 
R^ Rq, which allow to disregard the thermal noise in Polymer solutions at macroscopic 
scales, i.e. at scales much larger than typical interpolymer distances, can be described by a 
local elongation field R{x, t) which evolves according to 

^+u-VR=iVuf -R-^ (2) 

Taking the divergence of one easily sees that V ■ R decays in time and thus we can take 
R solenoidal. 

The effect of polymers on the fluid is in the modification of the stress tensor through 
an additional elastic component 11^ which takes into account the elastic forces of polymers 
Tlfj — {2vri / T){RiRj / R^) where v is the solvent viscosity and rj (proportional to polymer 
concentration) represents the zero-shear contribution of polymers to the total solution viscosity 
1^(1-1-77). The Navier-Stokes equation for the incompressible velocity field u{x, t) thus becomes 

du „ „ ^ 2vri R ■ VR 

— +u-Vu = -Vp+iyAu-\ — -2 — (3) 

at T Rq 

Equations Q and are the so-called uniaxial model for viscoelastic flows which can be 
also obtained starting from the linear Oldroyd-B model [11] by taking the limit of large 
elongations [12, 13]. Observe that by introducing the rescaled variable B = y^2^7//r(i?/_Ro), 
equations H2I3|I formally become the MHD equations for a plasma at zero resistivity with a 
linear damping —B/t [14,15]. In this representation the coefficient 77 disappears and thus 
the dynamics of the uniaxial model is independent of the concentration (which is physically 
consistent with the assumption of linearity and strong elongation). In the following we will 
thus take r] — 1 and absorb _Ro in the definition of i?: R ^ R/R^. 

The total energy (kinetic plus elastic) Et = {{u^) + (S^))/2 of the flow is dissipated at a 
rate 

-;.{|V X u\') - -{B^) - - Sr (4) 



dt 
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where is the viscous dissipation while the second term Sr represents the additional dissipa- 
tion due to the relaxation of polymers to their equilibrium configuration. 

In the following we will consider the statistics of stationary turbulent solutions of the 
viscoelastic model (|2I3|I integrated in a periodic box of size L — 2tt resolution 128"^ by 
means of a standard pseudo-spectral code for different values of the relaxation time t. For 
each r, a statistically stationary state is obtained by adding to ^ an external forcing term 
which acts on the largest scales by keeping their energy constant [16]. In stationary conditions 
the forcing injects energy with a mean rate £/ which balances the dissipation £/ = eu + £T- 
The turbulent regime for a viscoelastic fiow is controlled by two dimensionless parameters, 
the Reynolds number R\ = Urms^lv (where A = Urms/ {{dxUx)^)^^'^ is the Taylor niicroscale) 
and the Weissenberg number which is defined here as Wi — t/tk, where tk = {i^/SuY^^ is 
the Kolmogorov time. As a reference run, we integrated the standard Navier-Stokes equations 
((21 with 7] — 0, for which we have R\ ~ 87. In this Newtonian limit, we have also computed 
the Lagrangian Lyapunov exponent Xl which is a measure of the mean stretching rate. The 
dimensionless number A^t,, ~ 0.13 is consistent with known simulations [17]. The viscoelastic 
runs are performed for different values of the relaxation times correspoding to a Weissenberg 
number in the range 4.8 < W^i < 24 (i.e. 0.63 < Alt < 3.14). 
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Fig. 1 - Probability density functions of polymer elongations above the coil-stretch transition at 
Wi = 9.5 (squares), Wi = 12 (circles) and Wi = 24 (triangles). For Wi < Wi* = 6.5 the distribution 
is p{R) = 5{R). In the inset we show the growth of (R^) vs Wi. 

Although the uniaxial model is derived in the limit of strong elongation, it displays a clear 
coil-stretch transition as a function of Wi [18]. The limit Wi — > 0, i.e. t ^ 0, in |5J implies 
R = and Q recovers the usual Navier-Stokes equation for a Newtonian fluid. This coiled 
state persists until the mean stretching rate is comparable with the inverse relaxation time, 
Alt ~ 1 [19,20]. For larger values of Wi, the relaxation time of polymers is larger than the 
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smallest characteristic time of the tm-bulent flow and molecules start to be elongated. Figure^] 
shows the pdf of elongations p{R) at different values of Wi together with the mean square 
elongation. The coil-stretch transition at Wi* = 6.5 ± 0.5 (corresponding to Alt = 0.86) is 
evident. 

According to the Lumley criterion [18], as Wi increases above Wi* , polymers start to 
affect the dynamics of the turbulent cascade at the scale £l at which the eddy turnover time 
is of the same order of the relaxation time, i.e. £l ~ {er^Y^'^. For scales i > £l expect 
that the turbulent cascade is unaffected by polymers, i.e. with constant energy flux equal to 
the energy input e/, while for £ < £l we expect a reduced energy flux e — < £i because 
below the Lumley scale part of the flux is removed by elastic dissipation. 



0.25 




Fig. 2 - Newtonian (squares) and viscoelastic {Wi = 14.3, circles) spectra of kinetic energy, normalized 
with e^/'"^. The dotted line corresponds to Kolmogorov K41 scaling Exik) ^ k~^^^ . Lower inset: 
viscous dissipation normalized to the energy input ei as a function of Wi; the point size is of the 
order of the statistical uncertainty. Upper inset: drag coefficient D = siL/E^'^ versus Wi. 

By increasing Wi above Wi* two different scenarios are possible: the first is that elastic 
dissipation in |0J increases with Wi and energy flux at scales i < II vanishes (i.e. the Lumley 
scale £l becomes the new dissipative scale). A second alternative is that elastic dissipation 
removes only a fixed fraction of the flux. In this case Si, becomes independent of Wi and thus 
the energy cascade proceeds below (with a reduced flux). This latter scenario has been 
observed in shell models of viscoelastic fluids [21]. 

Our numerical simulations at increasing values of Wi indicate that the second scenario 
occurs. The inset of Fig. |2] shows that the ratio e,y/e/, which is by definition unity for 
Wi < Wi*, decreases for Wi > Wi* but already at Wi ~ 15 ~ 2Wi* saturates to a new 
value ~ 0.85. These results are in agreement with those reported in [5], where the reduction 
of vorticity in a viscoelastic solution was experimentally measured, and the Taylor microscale 
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A was observed to be practically unaffected by the presence of polymers. 

The above picture is supported by the comparison of the kinetic energy spectrum of a 
viscoelastic flow above the coil-stretch transition and the spectrum of the reference Newtonian 
flow. Fig.|21shows that, although the energy content at small scales is reduced by polymers, a 
power-law spectrum, characteristic of a turbulent cascade a la Kolmogorov, is present in the 
viscoelastic case as well. Moreover the effect of polymers on turbulence is local in scales, i.e. 
scales larger than are essentially not affected by the presence of polymers. Since the total 
kinetic energy Ek is dominated by large scales, we do not observe a significant variation of 
the "drag" coefficient, here defined a.s D — ejL/E^'^ (see the inset in Fig. |21 where it is shown 
that from Wi = to Wi = 24 we measure fluctuations of the drag of about 2% which are 
within the statistical uncertainty). This is at variance with the results reported in [7], where. 
However, it must be noticed that a different viscoelastic model was used, with a different large 
scale forcing mechanism. 
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Fig. 3 - Acceleration probability distributions for Wi — 0, 12, 24; the inner dotted curve shows a 
Gaussian for comparison. Inset: total arms, normalized with e^/^u'^^''', vs Wi. 



Since the turbulent cascade survives at scales smaller than ^l, we expect to observe some 
features of Newtonian turbulence in viscoelastic flows. One of the usual characteristics of 
small scale turbulence is the extremely intermittent acceleration which displays fluctuations 
much larger than the rms value [22,23]. In FigureOlwe report the statistics for the acceleration 
in viscoelastic turbulent flows at different values of elasticity above Wi* compared with the 
Newtonian case. In presence of polymers, the rms value arms is reduced with respect to the 
Newtonian case, again in agreement with experimental observations [24]. For sufficiently large 
values of Wi, arms becomes almost independent of Wi (and, at the largest value Wi — 24, is 
about 80% of the Newtonian case). This is consistent with a reduced value of energy flux at 
small scales shown in Fig. |21 indeed by compensating arms with the dimensional estimation 
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'^/^ this becomes virtually independent of Wi. 



The probability density functions shown in Fig. 13 indicate that relative fluctuations of 
turbulent acceleration are not affected by the presence of polymers. Indeed, the pdf of the 
rescaled quantity a /arms is found to be Wi- independent in all the range of Wi investigated. 
This is a remarkable result which, besides its intrinsic interest, has an important practical 
consequence: the fact that the shape of the acceleration pdf is not affected by polymers 
implies that it is possible to model small scale statistics in viscoelastic turbulence below the 
Lumley scale by means of the same models used for Newtonian fluids, by simply changing 
global quantities such as the energy flux. 

The acceleration a = du/dt in a viscoelastic flow has three different contributions from the 
rhs of (jSJ: pressure gradients, viscous and elastic contributions. In fully developed turbulence 
the viscous contribution is negligible and one may ask which is the dominant contribution for 
Wi > Wi*. Figure 21 (inset) shows that the pressure gradient contribution Up = ((Vp)^}^/^ 
(which is the only term in the Newtonian limit Wi = 0) is always dominant in the range of 
Wi investigated. Nevertheless, the contribution of polymers is not negligible: at the largest 
Wi — 24 the rms value of the elastic acceleration Uei = {2v/t){{R ■ Vi?)^)^/^ is almost 50% 
of the total acceleration arms ■ 
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Fig. 4 - Probability density functions of the cosine of the angle 6 between the vectors — Vp and 
R ■ Vii for Wi = 9.5 (squares), Wi = 14.3 (circles) and Wi — 24.0 (triangles). Inset: Leading 
contributions to the acceleration, normalized with the total rms acceleration, as a function of Wi: 
pressure gradient a^jarms (squares) and elastic stress contribution a^i/arms (circles). 



Since a ~ — Vp+ {2v/t)R ■ Vfi, the fact that ap ~ arms means that, increasing Wi^ the 
flow develops strong correlations between the pressure gradient and the elastic component in 
As a measure of this correlation we have computed the angle between the vectors — Vp 
and R ■ VR. Figure ^ displays the pdf of cos{9) and shows that, indeed, increasing Wi the 
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two vectors tend to be anticorrelated (i.e. cos{9) = —1) with higher and higher probabihty. 

In conclusion, we have studied the statistics of turbulent polymer solutions within the uni- 
axial model of viscoelastic flow. We have found numerical evidence for a coil-stretch transition 
at Wi* above which polymers affect the small scales of the turbulent flow. For Wi > Wi*, 
the energy flux is partially removed by polymer elasticity at the Lumley scale. This effect 
saturates and the turbulent cascade proceeds to scales smaller than the Lumley scale. As a 
consequence, small scale statistics, such as acceleration, display features typical of Newtonian 
turbulence. 
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